#include "dyn_quad_qGauss8.h"

qGauss8::qGauss8(Functor &_func, double _a, double _b)
    : qGauss(_func,_a,_b)
{
    populateVectors();
}

void qGauss8::populateVectors()
{
    weights.set_size(8);
    weights(0) = w4;
    weights(1) = w3;
    weights(2) = w2;
    weights(3) = w1;
    weights(4) = w1;
    weights(5) = w2;
    weights(6) = w3;
    weights(7) = w4;

    abscissas.set_size(8);
    abscissas(0) = -x4;
    abscissas(1) = -x3;
    abscissas(2) = -x2;
    abscissas(3) = -x1;
    abscissas(4) = x1;
    abscissas(5) = x2;
    abscissas(6) = x3;
    abscissas(7) = x4;
}
